Many-body interaction effects in doped and undoped graphene: 
Fermi liquid versus non-Fermi liquid 
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We consider theoretically the electron-electron interaction induced many-body effects in undoped 
('intrinsic') and doped ('extrinsic') 2D graphene layers. We find that (1) intrinsic graphene is a 
marginal Fermi liquid with the imaginary part of the self-energy, ImE(u), going as linear in energy 
uj for small u, implying that the quasiparticle spectral weight vanishes at the Dirac point as (lnw) -1 ; 
and, (2) extrinsic graphene is a well-defined Fermi liquid with ImE(oj) ~ tj 2 lnw near the Fermi 
surface similar to 2D carrier systems with parabolic energy dispersion. We provide analytical and 
numerical results for quasiparticle renormalization in graphene, concluding that all experimental 
graphene systems are ordinary 2D Fermi liquids since any doping automatically induces generic 
Fermi liquid behavior. 
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PACS numbers: 71.18.+y, 71.10.-w, 81.05.Uw, 73.63.Bd 

Recent progress in the experimental realization of a 
single layer of graphene [l[ has spawned tremendous in- 
terest and activity in studying the properties of this 
unique two-dimensional system. Graphene has a hon- 
eycomb real-space lattice structure, which comprises two 
interpenetrating triangular sublattices A and B, and a 
concomitant reciprocal space honeycomb structure with 
the hexagonal Brillouin zone cornering at the K points. 
At these so-called "Dirac points" the low-energy excita- 
tions satisfy the massless Dirac equation, and the band- 
structure for the "Dirac fermions" is uniquely differ- 
ent from the "Schrodingcr fermions" of the regular two- 
dimensional electron gas (2DEG) , exhibiting a linear en- 
ergy dispersion ej. = vk (having an effective "speed of 
light" v ~ 10 6 ms _1 [1]) with the conduction and valence 
bands connected at the Dirac point. 2D graphene is thus 
a rather unique chiral Dirac system. This apparently pe- 
culiar bandstructure of graphene is readily accounted for 
using a tight-binding model with nearest-neighbour hop- 
ping, where in the vicinity of the crossing of the energy 
dispersion relations for the two bands, the energy is lin- 
ear with respect to the momentum [2j. In an exciting 
recent development, this linear dispersion has been di- 
rectly observed using angle-resolved photoemission spec- 
troscopic (ARPES) measurements [3j. Another recent 
ARPES experiment Q finds that this linear spectrum 
manifests subtle many-body renormalization effects, pre- 
serving, however, the Landau Fermi liquid quasiparticle 
picture. On the theoretical side, there has been a num- 
ber of works in the literature on the effects of electron- 
electron interaction in graphene [H, [(|; however, the con- 
sideration has been solely based on undoped graphene so 
far, and explicit analytical results for either the doped or 
the undoped case are still lacking. The recent ARPES 
experiments are in doped graphene, and we show in this 
paper that doped and undoped graphene have qualita- 
tively different quasiparticle spectra. 

In this paper, we critically address the question of 



whether two-dimensional graphene is a Fermi liquid or 
not using the diagrammatic perturbation theory. We 
perform concrete theoretical calculations for the quasi- 
particle lifetime, renormalization factor, and effective ve- 
locity for both the doped and the undoped case. Wc 
find that extrinsic graphene manifests a Fermi liquid be- 
havior similar to the regular 2DEG, indicating that the 
quasiparticle description is indeed valid as substantiated 
by the recent experiment jij. For intrinsic (i.e. zero 
doping) graphene, we find that the quasiparticle lifetime 
scales linearly with energy oj above the Fermi energy 
(i.e. the Dirac point) and the renormalization factor 
at the Fermi energy vanishes, exhibiting a quintessential 
marginal Fermi liquid behavior. 

Thus, the quasiparticle behavior of 2D graphene de- 
pends crucially on whether the system is doped or not: 
While undoped intrinsic graphene is a marginal Fermi 
liquid, doped graphene with free carriers is invariably 
a garden-variety 2D Fermi liquid. Since the presence 
of charged impurities in the substrate (i.e. uninten- 
tional dopants) would invariably induce some carriers 
even in nominally undoped graphene, we conclude that 
the generic behavior of 2D graphene is likely to be that 
of a 2D Fermi liquid. 

The low-energy Hamiltonian for graphene is given by 
H = vcr ■ k, where v is the Fermi velocity of the Dirac 
fermions, <x is the set of Pauli matrices representing the 
two (A and B) sublattice degrees of freedom (through- 
out this paper we choose h = 1). This Hamiltonian 
describes a cone-like linear energy spectrum with con- 
duction band energy dispersion vk and valence band en- 
ergy dispersion —vk. The corresponding eigenstates are 
given respectively by the plane wave \k) multiplied by the 
spinors |±) = [1 ± e l ^]/\/2, where <f> = tan -1 (k y / l k x ) 
is the polar angle of the momentum k. The Green func- 
tion is given accordingly as Gk{ik n ) = Gk+{ik n )\+) (+| + 
Gh- {ik n )\— ){— |, where the label ±, also called the chi- 
rality, signifies the conduction band and valence band 



respectively, and Gk±(ik n ) = l/(ik n =p £fc) is the Green 
function in the diagonal basis, where we have denoted 
£fc = vk. The self-energy can also be similarly expressed 
as 

S fe (ife„) = E fc+ (ifc n )|+)(+| + E fc _(ifc n )|-)<-|, (1) 

where the leading-order self-energy in the diagonal basis 
is given by 



E fc ±(ifc„) = -k B T*^2 £ 

V q /e(q,iq n )] 



Gk+ q \{ik n + iq n ) 



1 ± Acos6> 



(2) 



with the Coulomb potential V q = 27re 2 /q, dielectric func- 
tion e and the scattering angle from k to k' = k + q 
denoted as 8. The factor (1 ± cos6>)/2 comes from over- 
lap of the eigenstates of the Hamiltonian H . Following 
standard procedure of analytic continuation, the retarded 
self-energy from Eq. ([2]) is obtained as 
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(3) 



where and ub are the Fermi and Bose distribution 
functions, respectively, and P stands for principal value 
integration. In the following, we proceed to calculate 
the quasiparticle lifetime, renormalization factor and ef- 
fective velocity at zero temperature for graphene in two 
distinct cases: extrinsic graphene with doping (Fermi en- 
ergy Bp > 0) and intrinsic graphene (ep — 0). Note 
that we calculate our self-energy (Fig. 1) in the leading- 
order infinite ring-diagram single-loop expansion approx- 
imation, which is essentially exact in the effective high- 
density regime of relevance to graphene. In fact, our self- 
energy calculation (Fig. 1) is a much better quantitative 
approximation to graphene than the so-called GW ap- 
proximation in regular metals and semiconductors which 
are never in the high-density regime. 

Extrinsic graphene. For extrinsic graphene, the quasi- 
particle is located in the vicinity of the Fermi level and 
we only need to consider the renormalized Fermi liquid 
parameters in the conduction band (or the valence band, 
depending on whether electrons or holes are the carriers 
as determined by doping). The lifetime of the quasi- 
particle — l/2r+ is obtained from the imaginary part of 
the self-energy Eq. ([3]) within the on-shell approximation 
£ = £k- 

Im£* + (e fe ) = \ £ £(1 + \cos8)V q [n B (\e k+q - e k ) 



A=± q 



-n F (\e k+q )\ Im 



e(g, Xe k+q - e k + i0+) 



(4) 
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FIG. 1: The Feynman diagrams for the screened Coulomb 
interaction within the RPA, and the retarded self-energy S R 
Eq. ([2}. The thin and thick wiggly lines denote respectively 
the bare and the RPA-screened interaction. The thin straight 
line stands for the bare Green function. This is the standard 
ring diagram approximation for self-energy exact in the r B <C 
1 limit. 



At zero temperature, it can be seen that interband scat- 
tering from the valence band A = — 1 vanishes in Eq. ([4|, 
and the expression reduces to (dropping the momentum 
k and energy e labels in the self-energy): 



Im£5 



p2ir pk 

-ire 2 / d(9(l + cos(9) / dk' k' 

JO Jk F 

Im[l/ £ (|fc-fc , |, gfc , -e k + iO+)] 
Vk 2 + fc' 2 - 2fcfc'cos(9 



(5) 



In order to maintain analytic tractability, we consider 
the long-wavelength x — qjlkp <C 1 limit and perform 
an analytical evaluation of the quasiparticle lifetime. For 
small x, the dominant contribution in Eq. ([5]) from the 
imaginary part of the dielectric function comes from low 
energies u — uj/vq <C 1, in which case the irreducible po- 
larizability for graphene (in the doped regime) is given 
in the leading order by Q II(q, w) ~ v{\ + iu) + (q/Av), 
where v — 4£^/27rt> 2 = 2kp/nv is the graphene den- 
sity of states at the Fermi level. We employ the ran- 
dom phase approximation (RPA) (Fig. [I} for the screened 
Coulomb potential and calculate the dielectric function 
as e(q,oj)/e m = 1 + V q Tl(q,uj) (here e m is the background 
dielectric constant). Substituting this expression for the 
dielectric function in Eq. (|5|), we obtain 



ImE^ = 

— / dk'k'{k 
Qtf Jk F 



(6) 
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Vfc 2 + k' 2 - 2kk'cos9 



where we have defined the Thomas-Fermi wavenumbcr 
as <7tf = Ae 2 kp/ve m . After carrying out the angular 
integral, expanding in the small parameter S = k — k' <C 
fcp, and integrating over k' ', we obtain 



P 2 

I mS R = _k. 

+ 87T£J 



\8epJ 2 



(7) 
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where = Sk — £f is the single-particle energy reck- 
oned from the Fermi level. Comparing with the corre- 
sponding expression for the regular 2DEG [8| Im£ R = 
(^/27T£ F )[ln(^ fe /16ei?) — 1/2], it is instructive to note 
that (1) the factor of two difference inside the logarithm 
comes from the linearity of the graphene spectrum; (2) 
the sign difference in the subleading term comes from the 
vanishing of the interband contribution from the valence 
band in Eq. {4J at zero temperature. 

Next we consider the renormalization factor for the 
Dirac quasiparticle Z — 1/(1 — A), where A is the deriva- 
tive of the real part of the self-energy Eq. (J2]) with respect 
to energy: 

13 f da; .., . . 

= ~2 9e / -^(jk+q\(ik n + IUJ) 

q,\ 00 

(l + \co S 9)V q /e(q,tu)\ k£=kFeF . (8) 

As usual, we perform the standard trick for evaluating 
the real part of the retarded self-energy by decompos- 
ing Eq. ([8]) into line and pole contributions [9J and per- 
forming analytic continuation ik n —> e + i0 + , and then 
integrating by parts the w-integral [l(J, obtaining 



dw Gk+ q \(iL<j)(l + \cos6)V q 



d 



1 



du> e(q, iui) 



(9) 
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The irreducible polarizability is n(q,w) ~ 1/(1 - 
u/^/u 2 — 1) + i(q/4v)/Vu 2 — 1. After some straightfor- 
ward algebra, we obtain the renormalization factor for 
the Dirac quasiparticle in the r s <C 1 limit as 



Z = l/ 



H 




•i)] 




7T V 





(10) 



where we have defined the interaction parameter r s by 
the ratio of the interparticle potential energy to the 
single-particle kinetic energy. It is interesting to note 
that whereas for regular 2DEG, r s = 2me 2 /fcffm is in- 
versely proportional to the square root of the 2D density, 
for graphene r s = e 2 /ve m is simply a constant ~ 0.73 
(with v ~ 10 6 ms _1 and e m ~ 3), indicating extrinsic 
graphene is essentially a weakly interacting (r s < 1) sys- 
tem. (We note in this context that ordinary metals have 
r s ~ 3 — 5 > 1, and the usual semiconductor-based 2DEG 
may have r s ~ 5 — 20 3> 1 !) In the leading order of r s , 
Eq. pop is similar to the corresponding expression for 
regular 2DEG [11] Z = 1 - (r s /27r)(l + tt/2). Beyond 
the leading order of r s , the q- integral in A has a logarith- 
mic divergence due to the interband contribution to the 
polarizability. Introducing a momentum cutoff k c of the 
order of inverse lattice spacing, we find the correspond- 
ing logarithmic correction ~ r 2 \a(k c /kp) appears in the 
second order of r s in the denominator of Eq. (fTUl) . 



For quasiparticles with regular quadratic spectrum, it 
is well-known Q that electron-electron interaction gives 
rise to effective mass renormalization of the Fermi liq- 
uid. In contrast, as the Dirac quasiparticle in graphene 
is massless, we have, instead, a renormalization of the 
quasiparticle velocity. We proceed to calculate the ef- 
fective velocity renormalization below, which is defined 
through v/v* = (1—A)/(1+B), where B is the derivative 
of the real part of the retarded self-energy with respect 
to momentum k and is obtained from Eq. (J8J) by replac- 
ing d/de — > (l/v)d/dk. Decomposing the expression for 
B into line and pole contributions and taking the ana- 
lytic continuation ik n — > e + i0 + , we find that the pole 
contribution vanishes and the expression for B becomes 



q.X 

(l + Xcos9)V q /e(q,iu)} 



I k,6 — kp ,6f 



(11) 



The evaluation of Eq. (fTTj) is similar to that for the renor- 
malization factor Eq. @, albeit more tedious. Carrying 
out the w-integral by integration by parts [lfj , and then 
performing the g-integral, we obtain the renormalized ve- 
locity at the Fermi level, within logarithmic accuracy, as 



v 
v 



1 



7T 



+ hi(r s ) 



4 kf 



(12) 



We also note that the same result here can be obtained 
more simply by using the static dielectric function and 
taking the derivative of the quasiparticle energy, consis- 
tent with the fact that the effective velocity at zero tem- 
perature in the lowest leading order in r s is only due 
to contribution from the static dielectric response. The 
first two terms in Eq. (I12p derive from the intraband 
contribution and are similar to the expression for the 
regular 2DEG [12] m*/m = 1 + (r s /2n)[2 + ln(r s /4)] 
whereas the last term arises solely from the interband 
contribution. In addition to v* at the Fermi level, we 
have evaluated v* at the Dirac point k = 0: v*/v = 
1 - r,{l + (l/4)ln[(l + 4r s )/4r s ] - (l/4)ln(fc c /fc F )}, high- 
lighting the fact that velocity renormalization is not uni- 
form for the entire spectrum but is in general a function 
of k. The renormalized spectrum therefore exhibits a 
small degree of nonlincarity imposed on the bare linear 
spectrum. 

Intrinsic graphene. The irreducible polarizability for 
intrinsic graphene is given by 0] 



IL(q,w) 



4 



0{vq 



6(uj — vq) 



\Jv 2 q 2 



v 2 q 2 



(13) 



Using Eq. fTB")) , the on-shell quasiparticle lifetime follows 
from Eq. Q as Im££ + (e k ) - £ q 9(\e k+q -e k \-e q )d(s k - 
£fc+q), which vanishes identically because of phase space 
restrictions imposed by the ^-functions. On the other 



hand, the imaginary part of the self-energy at k 
follows a linear relationship with u>, and is given by 



where 
f(r s ) 



lm££(fc = 0,w)=w/(r a ), 



[n(l-r s ) 



(14) 
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(nr s ) 2 7rr s - \J{irr s ) 2 - 4 



4^(tt7- s ) 2 - 4 ?rr s + ^(7rr s ) 2 - 4 



.(15) 



We emphasize that the linear relation is exact for all lo, 
indicating that intrinsic graphene is a marginal Fermi liq- 
uid. For arbitrary values of k we calculate ImH±(k,Lo) 
numerically, and Fig. 2 shows our calculated self-energy 
as a function of energy to for various values of momen- 
tum k. The numerical result for k — in Fig. 2 agrees 
precisely with the analytic result of Eq. (TT5|) . provid- 
ing a consistency check. As mentioned above we have 
ImE^: = for u> < et because there is no phase space 
available for virtual interband electron-hole excitations. 
However, for lo > et the virtual interband electron-hole 
excitations give rise to a finite ImX^ (i.e. finite quasi- 
particle lifetime). Initially ImE^ (Im£_) rises sharply 
(slowly), and for large values of momentum they increase 
linearly with the same slope as that for ImE^/e = Q,lo) 
(c.f. Eq. (fT5|) ). Note that there is no plasmon con- 
tribution to the imaginary part of the self-energy for 
intrinsic graphene. The contribution of the interband 
electron-hole excitations gives rise to the linear behavior 
of ImE^:. However, for doped graphene the contribu- 
tions of the interband electron-hole excitations is com- 
pletely suppressed due to phase space restrictions at zero 
temperature. The contributions of the intraband virtual 
single-particle excitations and/or the virtual excitations 
of plasmons give rise to higher powers of lo (i.e. lo 2 ) in 
the imaginary part of the self-energy in the doped case, 
which restores the usual Fermi liquid behaviour. Thus, 
the qualitative difference between intrinsic (ImE ~ lo)) 
and extrinsic (ImE ~ lo 2 ) graphene can be completely 
understood by noting that the intrinsic system is an 
insulator (albeit a zero-gap semiconductor with no in- 
traband single-particle excitation) and the extrinsic case 
has a Fermi surface with intraband single-particle excita- 
tions. Thus, any doping of graphene (intentional or unin- 
tentional) will immediately suppress its marginal Fermi 
liquid intrinsic character, converting it to a regular 2D 
Fermi liquid. 

The renormalization factor can be obtained from 
Eq. ([9]), whereupon evaluating we find that A ~ J dq/q 
diverges logarithmically, which is due to the unscreened 
nature of the Coulomb potential (the same divergence oc- 
curs for the exchange energy of a regular 2DEG). There- 
fore the renormalization factor Z = 0, showing that 
as doping goes to zero, the magnitude of the step at 
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FIG. 2: (Color online) Imaginary part of the self-energy for 
the intrinsic graphene, Im£+(fc, lo) (solid lines), ImE^(/;,a)) 



(dashed lines), for different 
0, 1, 2 eV/hv (black line '0', 
spectively). Note that for k = 



values of momentum, k = 
red line '1', blue line '2', re- 
0, ImE5 = ImE H . 



the Fermi energy e = also shrinks to zero, approach- 
ing the Dirac point where the notion of a "Fermi sur- 
face" no longer applies, the quintessential behaviour of 
a marginal Fermi liquid. Using the Kramers-Kronig re- 
lations and from the fact that Im£+(fc = 0,lo) ~ ui, 
we have ReE^O = 0, lo) ~ wlnw, and the renormal- 
ization factor Z 1/lncj, which approaches zero log- 
arithmically as lo — > at the Fermi energy. More- 
over, the spectral function p{k — Q,lo) = ImS^(fc = 
0,w)/{[ImS^(fc = 0,lo)} 2 + [u- ReSf(k = 0,to)} 2 } di- 
verges as p(k — 0,lo) r~j 1/lo (lnw) 2 . In addition, we also 
find that the effective velocity v* /v = l + (r*/4)ln(fc c //c F ) 
(here r* = r s /(l + Ttr s /2) is the renormalized interaction 
parameter) diverges as doping kp — ► 0. These results for 
the intrinsic graphene are consistent with Ref. 0], where 
the renormalization group approach is used to arrive at a 
similar conclusion. We note in passing that the case for 
a purely undoped 3D system with agapless linear energy 
dispersion was considered in Ref. 1 1 31 ] . and was found to 
exhibit marginal Fermi liquid behavior with a logarith- 
mic energy dependence in ReS comparable to our results 
for intrinsic graphene. 

In conclusion, we have presented a calculation, for- 
mally exact in the r s -C 1 limit, for the renormalized 
Fermi liquid parameters for both extrinsic and intrinsic 
graphene. We find that for extrinsic graphene the ana- 
lytical results for the quasiparticle lifetime, renormaliza- 
tion factor and effective velocity show no deviation from 
the usual Fermi liquid behaviour, and the Fermi liquid 
description is robust. On the other hand, with precise 



5 



zero doping, intrinsic graphene exhibits a quasiparticle 
lifetime linear in the excitation energy and a zero renor- 
malization factor, indicating that the Fermi liquid de- 
scription is marginal at the Dirac point. With a finite 
Fermi energy in the extrinsic graphene, the interband 
single-particle excitations which give rise to the linear in- 
dependence of the quasiparticle lifetime (and hence the 
marginal Fermi liquid behavior) in the intrinsic graphene 
are suppressed, bringing the system back to a usual Fermi 
liquid. Since some finite doping is invariable in real sys- 
tems, we predict real 2D graphene to be generically a 
Fermi liquid. 
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